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Abstract 



^: 

, We investigate the effect of different mass of a Bose- and a Fermi-particle on the collective 
oscillations of the degenerate boson-fermion mixtures. In particular we consider the monopole and 
t.he quadrupole modes of the oscillations and study their characters and the frequencies by using 

5-H ; 

^ variational-sum-rule approach. We find that for both the modes there exists a critical value of 
O the ratio of boson-fermion mass below and above which the character and the frequency of the 
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espective modes are significantly different. 
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I. INTRODUCTION 



In recent years many theoretical and experimental studies involving dilute mixture of 
trapped ultra cold gases of bosonic and fermionic atoms have been reported in the literature. 
These mixtures provide a convenient way to achieve degenerate fermionic gas by means 
of sympathetic cooling as the conventional evaporative cooling methods used to obtain 
Bose-Einstein condensates (BEC) are not applicable to fermions. Besides this Bose-Fermi 
mixtures can be used to study many aspects of quantum statistics. Experimentally, stable 
3EC immersed in a degenerate Fermi gas have been realized with '^Li in ^Li ll, '^^Na in ^Li 
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2|, and ^^Rb in "^^K j3|, y? ]. The first two Bose-Fermi mixtures are characterized by positive 
inter-species scattering lengths or repulsive interaction between bosons and fermions. On the 
other hand, mixture of ^'^Rb and '^^K give rise to large attractive boson-fermion interaction. 
Recently, by exploiting Fesbach resonance the inter-species interaction has been tuned with 
a large range positive and negative values. 

Theoretically, the trapped system of Bose-Fermi mixture have been studied within the 
mean-field approximation to determine the bosons and fermion density profiles at zero tem- 
perature. These studies mainly focused on the phenomena of phase separation for repulsive 
and collapse mixture for attractive boson-fermion interaction respectively have been studied 
extensively [g], Q, Is], y, 10, [ll, 12, Q, Q. Beside ground state density profiles the dynamic 
properties of Bose-Fermi mixture like free expansion [15l and the spectrum of collective ex- 
citations Q, 0, Q, Q, Q have also been studied. All these studies have clearly showed 
that the spectra of collective excitations bear unambiguous signatures of phase transitions 
that Bose-Fermi mixtures undergo under appropriate conditions. Moreover, for BEC it is 
well established that the frequencies of collective oscillations can be measured with very 
high precision. Consequently, it is natural to expect that by measuring the frequencies of 
the collective oscillations of a trapped Bose-Fermi mixture it would be possible to determine 
the parameters at which mixture undergo phase separation or collapse. 

Keeping this in mind we focus our attention on the effect of unequal mass of a Bose- 
(mfe) and a Fermi-particle {mf) on the frequencies of collective oscillations. At this point 
we note that most of the studies on the collective oscillations assumed rrih = TTif. In reality, 
for all the experimental realizations as mentioned above, m^, ^ rrif. For example, in the 
case of extensively studied Rb-K mixture, mb/mf ^ 87/40. In this paper we show that 
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rrih ^ rrif leads to siErnificant alteration of the nature and the frequencies of the modes of the 
collective oscillations over the equal mass case [16|, |19|. In order to calculate the frequencies 
of collective oscillations we adopt well established sum rule approach of many-body response 



theory 
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22j. It has already been shown that the sum rule approach yields quite accurate 



results for the frequencies of collective oscillations of trapped Bose-Fermi mixtures [17 |. 
Using sum rule approach, we derive analytical expressions for the frequencies of monopole 
and quadrupole modes of oscillations of a Bose-Fermi mixture with mass of a Bose-particle 
being different from that of a Fermi-particle. Our expressions are generalizations of the 
results of Ref . [l^ . These analytical expressions for frequencies are then used to determine 
critical mass ratio (m{,/mj)c) at which character of the modes undergoes a change resulting 
in the modifications of the frequencies significantly. 

The content of the paper is as follows. In the next section we present the derivations of 
analytical expressions for the frequencies followed by a discussion of our results in Section 
III. The paper is concluded in Section IV. 



II. VARIATIONAL SUM RULE APPROACH 



The basic result of the sum-rule approach [2l|, |22| is that the upper bound of the lowest 
excitation energy is given by 

Mlex = a/—, (1) 
V mi 

where 

mp = Y,mF\n)\^hjnor, (2) 

n 

is the p-th order moment of the excitation energy associated with the excitation operator 
F and Qex is the frequency of excitation. Here huno = En — Eq is the excitation energy of 
eigenstate \n) of the Hamiltonian H . Moreover, Eq. ([1]) can be used for computation of the 
excitation energies by exploiting the fact that the moments mi and can be expressed 
as expectation values of the commutators between F and H in the ground state |0). These 
relations are 

mi = i(0| [F\[H,F]] |0), 

= \{^\[[F\H],[[H,[H,F]]]\Q). (3) 
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The main advantage of the sum-rule approach is that it allows us to calculate the dynamic 
properties like excitation frequencies of many-body systems with the knowledge of ground 
state |0) (or the ground-state density) only. 

For the purpose of calculation one needs to choose an appropriate excitation operator 
F. In this paper we consider the monopole and the quadrupole modes of the collective 
oscillations. The excitation operators corresponding to these modes are defined by 

^.^ = /^(r)±/i(r), (4) 

where a stands for monopole or quadrupole modes, the indices h/ f denote bosons/fermions 
and the functions i^are given by /„ = to. monop.e and = 3.^ - <or ,„ad.npole 
modes. Following [16|] we take the excitation operators as a linear combination of and 
F~ to simulate the effect of mixing of boson and fermion oscillations. It is written as 

Fa{9) = F+cosd + F-smO (5) 

where the mixing angle 9 lies between — 7r/2 and 7r/2. The character of the mode is given 
by the value of 9, for example, = 7r/4 for the bosonic and —ti/A for the fermionic modes, 
6* = for the in-phase oscillation, and 9 = -njl for the out-of-phase oscillations of the two 
types of particles. The value of mixing angle for each mode is determined by minimizing the 
corresponding frequency. The presence of two kinds of particle leads to two types of collective 
oscillations for each multipole which are orthogonal to each other. Out of these two modes 
the low-lying mode is excited by the operator of Eq. ^ and the operator F^sin9i —F^cos9i 
for each a excites the orthogonal high-lying mode. To study the effect of unequal mass we 
focus our attention on the low-lying mode only. 

In order to calculate the moments mi and m^, we use the expression for the energy 
functional (in terms of boson pb and fermion pf densities) of a harmonically trapped Bose- 
Fermi mixture as given by (23I 

E[pb,pf] = TM + VM+E':'Jpb] 

+ Tf[pf] + Vf[pf]+Eli[pb,pf] (6) 

where Tb[pb] (H[Pfe]) and Tf[pf] (V/[p/]) represent the kinetic energy (the trapping energy) 
functionals of bosons and fermions respectively. The remaining two terms E'^^^lpb] and 
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^intiPbj Pf] represent the boson-boson and the boson-fermion interaction energy functionals 
respectively. In writing above energy functional we have neglected the fermion-fermion 
interaction energy term as fermions are assumed to be spin polarized. The kinetic energy 
functional T(,[pb] is given by 

TM = 1^1 l^VP'bld'r. (7) 
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On the other hand, exact form of the kinetic energy functional for fermions is not known. In 
this paper we work within the Thomas-Fermi approximation and use following expression 
for the kinetic energy functional 

TM' — i^-^rJrrMd'r. (8) 

This approximation is quite accurate when the number of fermions is fairly large and in this 
paper we choose number of fermions (Nf) such that the TF approximation remains valid. 
The harmonic trapping energy functional for spherically symmetric trapping potential are 
given by 

Vf[Pf] = Imfuj j r'pf{r)£r, (9) 

where ujf, and Uf are frequencies of the trapping potential for bosons and fermions respec- 
tively. In accordance with the experimental settings we use cof — ■<^mb/mfUJb for the purpose 
of our calculations. Finally the forms for the interaction energies within the mean-field ap- 
proximation are given by 

E^ntiPb] = 9bb J Pl{v)dh 

E\L[pb.Pi] = 9bf j Pb{v)pf{v)dh (10) 

The boson-boson coupling strength is given by g^i, = Airfi^aijb/mh, where 055 is the 
boson-boson s-wave scattering length. The boson-fermion coupling strength reads Qbf = 
Airh'^abf/mr, where Ubf is the boson-fermion s-wave scattering length and = mbruf / {mb-\- 
TTif) is the reduced boson-fermion mass. 

Before proceeding further we wish to note here that representation of ground-state en- 
ergy as functional of density is ensured by Hohenberg-Kohn theorem of density functional 



theory [2J]. The ground-state densities ph and pf can be determined by imposing stationary 
conditions: 

^E[pb,pf] _ _ 5E[pi,,pf] 



11) 



Spb Spf 
where pf, and pf are boson and fermion chemical potentials, respectively. The chemical 

potential are fixed by the normalization conditions / pb{r)d'^r = Ni, and / pf{r)d^r = Nj, 

where Nh and Nf representing number of bosons and fermions, respectively. In this paper 

we exploit this variational nature of the ground-state energy to determine the ground-state 

density profiles of bosons and fermions. More about this approach will be discussed later in 

this section. 

Using the above energy functional (Eq. ([6])) along with the Eq. ([3]) we find after a tedious 
although straightforward algebra following expression for the frequency of the monopole 
mode of collective oscillations 



UJ, 



{Xm + + + 2 siu 9 COS 9 {X^ — ^m) — 2 siu^ 9Z„ 
{Am + Bm) + 2 sin 9cos9 {A^ - Bm) 



'12) 



with 



Xrr, 

Y 

Bm 



m + m 

, 2 

rrif 



9Et + 9EZ + ^Af-A 



int 



int 



ATf + AVf + + 3At, - A 



bf 



int 



2 



Vf 



(13) 



In the above expressions Tbj, Hj, E\^^ and E\1^ are the respective integrals defined in Eqs. 
© - (|TOl) evaluated with the ground-state boson and fermion densities. Besides these energy 
integrals, remaining three quantitities A/, A?,, and A are given by 



Afe 

A/ 
A 



9bf 
9bf 
9bf 



rpbir)'-PMdh 
dr 

rpf{r)^^d\ 
dr 

dpb{r) dpf{r] 



d^r. 



dr dr 

The expression for the frequency of the quadrupole mode is given by 



{Xg + Yg + Zq) + 2sin9cos9 (X, - Yq) - 2sin^9Zg 



{Aq + Bq) + 2sin9cos9 {Aq 
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(14) 



(15) 



with 

and Aq = A^, Bq = Bm- 

The expressions given by Eqs. (fT2|) - (ITSl) constitute the main results of this paper and 
they correctly reduce to the corresponding equations of Ref. 1^ for the case of mb/mf = 1. 
In order to calculate the frequencies, we first need to determine the ground-state densities 
of bosons and fermions. As mentioned before, we accomplish this task by employing a 
variational approach similar to the one described in Ref. |8|. For fermions we apply TF 
approximation and within this approximation the density of fermions is given by 

PA^)=[^y'-^2if^f-Vf{r)-g,M-)f"- (17) 

On the other hand, the density of bosons is determined by choosing appropriate variational 
forms for it. In particular for negative values of boson-fermion scattering length abj, we use 
single-parameter gaussian form for the boson density given by P6(r) = A;,e~"^^, where A^ is 
the normalization constant fixed by the normalization condition J pb{r)d^r = Ni, and a is 
the single variational parameter determined by minimization of the energy functional of Eq. 



([6]). It has already been shown 25j that gaussian ansatz leads to incorrect density profile 
for positive values of a^j specially when atf ^ CLbb- In order to circumvent this problem, we 
employ a two-parameter form 



n 

proposed by Fetter [26] to study BEG. In Eq. (ITSl) . 9{R — r) is the step function, R and 
A are the two variational parameters again determined by minimization of the total energy. 
Here we mention that the ground-state calculations performed with this variational approach 
yields results which are quite close to the corresponding results reported in Refs. [q. 16| 

By using the ground-state densities obtained by variational approach we evaluate the 
required integrals to calculate the frequencies of the collective oscillations given by Eqs. 
( fT2|) and (|T5|l . The results of these calculations are discussed in the next section. 



III. RESULTS AND DISCUSSION 



We begin this section with the discussion of results for the monopole mode of the collective 
oscillations followed by the results for quadrupole mode. Before discussing the results we note 
that in this paper all the numerical calculations are performed for Nf, = 10^ and Nf = 10^. 
Following Ref. [l6| the boson-boson scattering length is chosen to be abb = CLbb/h = 0.2/ (Svr) 



( where lb = ^yh/nibUJb ) and the boson-fermion scattering length is varied from a negative 
value characterizing collapse of the mixture to a large positive value at which Bose-Fermi 
mixture undergoes phase separation. Such wide variation of the boson-fermion interaction 
strength has recently been achieved by tuning the magnetic field at a Fesbach resonance [5] 
We wish to emphasize here that the above choice of the parameters does not lead to any 
loss of generality of the results presented in this paper. 



A. Monopole mode 

In order to demonstrate the effect of different mass of Bose- and Fermi- particles, we first 
display the frequency and the mixing angle of the monopole mode as a function of coupling 
strength k = abj/abb in Fig. 1 for two different values of mass ratio, namely mb/nif = 1 
and nib/mf = 87/40. The latter ratio corresponds to the Bose-Fermi mixture of ^'^Rb and 
'^^K atoms. Fig. 1 clearly shows that the monopole mode of the collective oscillations 
exhibits significantly different behaviour when the ratio rrib/mf is changed from 1 to 87/40. 
The differences in the behaviour of the mixing angle which characterizes the nature of the 
mode of the collective oscillations are clearly visible even in the weak interaction regime 
(k < 1). For example, in the weak boson-fermion interaction regime the monopole mode 
is dominantly fermionic in nature for nib/mf = 1, whereas it becomes dominantly bosonic 
for mb/rrif = 87/40. Consequently, the frequency of the monopole mode shows completely 
different behaviour as the mass ratio is changed. This result has motivated us to investigate 
the effect of different mass of a Bose- and a Fermi-particle on the collective oscillations of 
boson-fermion mixtures. At this stage it is natural to ask following question: is it that 
any mass ratio different from unity will change the behaviour of monopole mode of the 
collective oscillations or there exists a critical value of mass ratio {mb/mf)c at which change 
in the behaviour of monopole mode occurs? In order to answer this question, we focus our 
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attention on the collective oscillations of non- interacting {abf = 0) boson- fermion mixture 
as we have already seen that effect of non-unity mass ratio is appreciable even for very small 
values of k around zero. Moreover, for non-interacting Bose-Fermi mixture one can derive 
analytical results for the frequencies and the mixing angle, which makes the analysis more 
lucid. By using the expression for the frequency of the monopole mode given by Eq. (1121) for 
non-interacting case the condition of extremum dujm/dO = leads to two possible solutions 
for 6: 6 = ±7r/4. Now to determine the value 6min out of the above two values at which 
becomes minimum, we impose the condition {(Pujm/ dO"^) < for 9min- This condition then 
leads to 



Gmnr:. slu 26'™,i„. < 0, whcrc G 



(19) 



From the above inequality, we infer that when Gmon is positive (negative), the correspond- 
ing mixing angle of the monopole mode becomes 6min = — 7r/4(7r/4). In the limit of TF 
approximation for bosons (Tb/Vb = 0), Eq. (ITUi) reduces to 



G'mon sin 2emin < 0, where 



5-4 



ruf 



(20) 



For rub = TUf^^S considered in the previous works on the collective oscillations of Bose-Fermi 



mixtures 
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19| , Gmo^re Eq. fl2Ql) is positive and consequently the mixing angle minimizing 
the monopole mode is —it /A. Therefore, the monopole mode for a non- interacting Bose-Fermi 
mixture with nib = rrif is purely fermionic in character and the frequency of the oscillation 
is given by = ^Ub- These results are clearly illustrated in Fig. 1 by the dashed curves. 
On the other hand, for rUb 7^ nif, it can be easily inferred from Eq. ( l20l) that Omin can 
assume value —tt/A or 7r/4 depending on the value of the ratio mb/rrif. The transition from 
^min = — 7r/4 to 9min = 7r/4 will take place at a critical value {mb/mf)c = 5/4 = 1.25. For 
mb/rrif < 1.25, G^ara is positive and consequently 9min = — vr/4 indicating that the monopole 
mode is purely fermionic in nature similar to the case of rrib = rUf . In contrast to this, for 
mb/rrif > 1.25, the functin G^^ is negative and the Eq. (1201) is satisfied for 9jnin = 7r/4 
leading to a purely bosonic monopole mode with the frequency cUm = V^ujb- The results 
shown in Fig. 1 with solid curves correspond to this situation as mb/mf = 87/40 > 1.25. 
When Tb/Vb 7^ the critical mass ratio is given by 

'mb\ 5 Tb 



mjj^ 4 4H' ^^^^ 



and to determine the value of 6min for general case, we plot Gmon as a function of the mass 
ratio mb/mf in Fig. 2. It can be seen from Fig. 2 that even for the general case the value 
of the critical mass ratio is very close to TF case result (mf,/mj)c = 1.25. Therefore, from 
the above results we infer that depending on the value of mass ratio rrib/mf the monopole 
mode of the boson-fermion mixture can be dominantly bosonic or fermionic in character. 
The transition from one to other occurs at critical mass ratio (mft/m/)c = 1.25. 

Having found the nature of the monopole mode for the non-interacting Bose-Fermi mix- 
ture, we next focus our attention on the interacting case {a^f ^ 0). For the interacting 
case, we obtain the value of 6min by numerically minimizing the monopole mode frequency 
Urn- The results of these numerical calculations are also displayed in Fig. 1. Wc find that 
for large values of boson-fermion interaction strength both in the positive and the nega- 
tive directions the monopole mode becomes coherent superposition of bosonic and fermionic 
oscillations. The proportion of mixing of these two types oscillations depends on the boson- 
fermion mass ratio and the interaction strength k. The monopole modes corresponding to 
both mb/mf < {mb/mf)c and mb/mf > {mb/mf)c cases tend toward in- phase oscillations 
{d,nin = 0) of the Bose-Fermi mixture as it approaches either the collapse or the spatial phase 
separation regimes. Wc note here that depending on the ratio of the Bose-Fermi mass the 
monopole mode follows either the lower branch characterized by purely fermionic oscillations 
at zero interaction to the coherent superposition of bosonic and fermionic oscillations for 
finite values of interaction strength or the upper branch representing purely bosonic to the 
coherent superposition of bosonic and fermionic oscillations. Consequently, the frequency of 
the monopole mode show drastically different behaviours for mass ratios below and above 
the critical limit. Moreover, from the behaviour of the monopole mode discussed above, we 
note that for a given mass ratio the character of the monopole mode at zero boson-fermion 
interaction strength can be employed to predict the nature of this mode at finite values of 
the interaction strength. For example, another experimentally achieved Bose-Fermi mixture 
composed of '^Li and ^Li has mass ratio mb/mf — 7/6, which is less than 1.25. From the 
above discussion, we infer that at zero boson-fermion interaction strength monopole mode 
is purely fermionic in character and for finite values of interaction strength it will follow the 
lower branch similar to the case of mb/mf = 1 as shown by dashed lines in Fig. 1. 
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B. Quadrupole mode 



Now we consider the quadrupole mode of the collective oscillations of a trapped Bose- 
Fermi mixture and carry out an analysis similar to the one discussed above for the monopole 
mode. Using Eqs. (ITB]) and (ITBi) for dhj = 0, we find that for quadrupole mode also 9 = — 7r/4 
and 6 = 7i/A make the frequency Uq extremum. The condition for obtaining mixing angle 
Omin that minimizes the frequency of the quadrupole mode is given by 



Gguad sin 26'™„ < where Gguad 
Within the TF approximation for the bosons above criterion reduces to 



(22) 



Gquad^^T^'^^min < whcrC G^f^^ 



1-2 



sin29mAr,.. (23) 



For mfe = mj the function G^^^ in Eq. (123!) is negative and the minimum condition is 
satisfied if 6min = vr/4. Therefore, for non-interacting Bose- Fermi mixture the quadrupole 
mode is of purely bosonic character and the frequency of the oscillations is cUg = V2uJb [16]. 
However, for ^ mj, the quadrupole frequency can become minimum either for Omin = 
7r/4 or 9min = —vr/4 depending on whether mb/rrif > {mb/mf)c or nib/mf < {mb/mf)c 
respectively, where {mb/mf)c = 0.5. Thus, for rrib/mf < 0.5, the quadrupole mode shows 
purely fermionic character and the corresponding frequency of the non-interacting Bose- 
Fermi mixture is given by Uq = 2ujb. On the other hand, when mb/mf > 0.5, the quadrupole 
mode possesses purely bosonic character with frequency Uq = V2uJb- 

For the general case {Tb/Vb ^ 0) the critical mass ratio at which quadrupole undergoes a 
change in character is given by 

In order to determine the mixing angle that minimizes the quadrupole frequency for the 
general case, we plot in Fig. 3 the Gquad as a function of mb/mf. Again it can be seen that 
corresponding to A*";, = 10® the critical mass ratio for quadrupole mode is very close the TF 
result {mb/mf)c = 0.5. 

Like monopole mode, the mixing angle and the frequency of the quadrupole mode at finite 
value of the boson-fermion interaction strength are obtained numerically by minimizing the 
quadrupole frequency. The results of such calculation for two different values of mass ratio 
namely, (mb/mf) = 1/3 and (mb/mf) = 1 are shown in Fig. 4. These two figures clearly elu- 
cidate the results on the quadrupole mode discussed above. For example, for (mb/mf) = 1/3 
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the minimum criterion Eq. (l22|) is satisfied by 9min = — vr/4 corresponding to a fermionic 
quadrupole mode with the frequency Uq = 2ujb at zero value of the boson-fermion interaction. 
For finite values of the boson-fermion interaction strength the quadrupole mode no longer 
remains a purely bosonic or fermionic in character, rather it becomes a coherent superposi- 
tion of both. The proportion of the mixing two modes depends on the boson-fermion mass 
ratio and the interaction strength. Like monopole mode the superposed quadrupole mode 
either lie on the lower or the upper branch depending on the mass ratio at finite values of 
the boson-fermion interaction strength. As a result of this the frequency of the quadrupole 
mode show different behaviour for mb/mf < {mb/mf)c and mb/rrif > {mb/mf)c and this is 
clearly elucidated in Fig. 4. 



IV. CONCLUSION 



In this paper we have studied the effect of different mass of a Bose- and a Fermi-particle 
on collective oscillations of trapped Bose-Fermi mixtures. In particular we have considered 
the monopole and the quadrupole modes of collective oscillations and investigated the effect 
of unequal mass on the characters and the frequencies of the modes. The calculations of 
frequencies are performed by employing sum-rule approach of many-body response theory 
in conjunction with a variational method to obtain the ground-state densities of bosons 
and fermions. By considering the non-interacting Bose-Fermi mixtures we have determined 
the critical mass ratio both for the monopole and the quadrupole modes. At these critical 
values of mb/mf the monopole and quadrupole modes undergo change in character and 
consequently the frequencies. In order to study the effect of unequal mass for the finite values 
of the boson-fermion interaction strength we determine the mixing angle characterizing 
the nature of the mode by numerical minimization of the corresponding frequencies. We 
have shown that the characters of both monopole the quadrupole modes of a Bose-Fermi 
mixture are decided by their corresponding nature at non-interacting regimes. The recent 
achievement of trapped degenerate Bose-Fermi mixtures of atomic gases in laboratory and 
the capability of tuning the boson-fermion interaction strength by Fesbach resonance make 
the experimental verification of the results presented in this paper feasible. 
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Figure captions 

Fig.l The frequency (upper part) and the mixing angle (lower part) of the monopole 
mode as a function of the boson-fermion interaction strength k for two different values of 
mh/nif : mf,/mf = 1 ( dashed curve) and mh/nif = 87/40 ( solid curve). The horizontal 
upper and the lower dotted lines in the bottom panel denote the mixing angles for pure 
bosonic and fermionic modes respectively. 

Fig.2Plot of function Gmon of Eq. (fT9|) as function of the ratio nih/mf. The intersection 
with the dotted horizontal line denotes the critical value of ratio of masses of a boson and 
a fermion (mf,/mj)^. 

Fig.SPlot of function Gquad of Eq. fl2^ as function of the ratio mfy/mf. The intersection 
with the dotted horizontal line denotes the critical value of ratio of masses of a boson and 
a fermion {rrih/mf)^. 

Fig. 4 The frequency (upper part) and the mixing angle (lower part) of the quadrupole 
mode as a function of the boson-fermion interaction strength k for two different values of 
rrib/mf: rrib/mf = 1 ( dashed curve) and mb/nif = 1/3 ( solid curve). The horizontal upper 
and the lower dotted lines in the bottom panel denote the mixing angles for pure bosonic 
and fermionic modes respectively. 



15 



